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Abstract 

We use the Raychaudhuri equation to probe certain aspects related to the gravitational 
collapse of a charged medium. The aim is to identify the stresses the Maxwell field exerts 
on the fluid and discuss their potential implications. Particular attention is given to those 
stresses that resist contraction. After looking at the general case, we consider the two oppo- 
site limits of poor and high electrical conductivity. In the former there are electric fields but 
no currents, while in the latter the situation is reversed. When the conductivity is low, we 
find that the main agents acting against the collapse are the Coulomb forces triggered by the 
presence of an excess charge. At the ideal Magnetohydrodynamic (MHD) limit, on the other 
hand, the strongest resistance seems to come from the tension of the magnetic forcelines. In 
either case, we discuss whether and how the aforementioned resisting stresses may halt the 
contraction and provide a set of conditions making this likely to happen. 

PACS numbers: 04.20-q, 04.40.Nr, 95.30.Qd 

1 Introduction 

Gravitational collapse has been one of the main topics of research interest in relativistic astro- 
physics for many decades. The implications of nonzero charges, together with those of electro- 
magnetism, for the ultimate fate of a self-gravitating massive body, also have a long research 
history. Over the last fifty years, or so, several authors have looked into this question from a va- 
riety of perspectives and by means of different techniques. Some of the available work addresses 
static problems, while other authors have considered dynamical situations. For a representative 
list the reader is referred to [T]-|21j. with more references given in [19J. A number of these stud- 
ies have raised the possibility that a contracting charged medium may 'bounce back' instead of 
'falling' into a singularity. The first signs were probably those seen in Melvin's 'magnetic uni- 
verse' and in the studies of collapsing charged dust by Novikov [2]-[6]. With very few exceptions, 
all of the aforementioned work uses metric-based techniques. 

An alternative approach to gravitational collapse is to employ covariant methods and thus in- 
volve the Raychaudhuri equation [22j. The latter determines the average contraction/expansion 
of a self-gravitating medium, by monitoring the cross-sectional area of the particle worldlines. 
Raychaudhuri's formula is the key equation of gravitational collapse and has played a central 
role in the development of the singularity theorems of the late 1960s and early 1970s (for recent 
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reviews and extensive discussion the reader is referred to |23|-|26|). A major difference in the 
nature of the Raychaudhuri equation, when applied to a charged environment, is that the world- 
lines of the matter are no longer geodesies. This happens every time non-gravitational forces are 
included into the system and the particles' motion is not dictated by gravity alone. In our case, 
the fluid kinematics are determined by the combined action of the Einstein and the Maxwell 
fields. As a result, new stresses are added to Raychaudhuri's formula. These are primarily due 
to the Lorentz force and depend on the electric properties of the charged medium. On whether, 
in particular, the conductivity is low or high. In the first case, there are no electric currents 
despite the presence of a finite electric field. At the ideal Magnetohydrodynamic (MHD) limit, 
on the other hand, the situation is reversed. 

For our purposes the stresses of interest are those that tend to resist contraction. When 
dealing with poorly conductive media, we find that the main resistance comes from the familiar 
Coulomb forces triggered by the presence of a net charge. The overall effect, which is partly 
Newtonian and partly relativistic, depends on the specific charge of the matter. When the latter 
exceeds a particular threshold, the Coulomb stresses become strong enough to outbalance the 
gravitational pull of the matter fields and bring the collapse to a halt (see § 14.11 14. 2D . Analogous 
dynamical results, both qualitatively and quantitatively, have also been obtained through semi- 
covariant and metric-based methods [U [16] . When the conductivity is high the electric fields 
vanish and the currents keep the magnetic forcelines frozen-in with the fluid. Then, it appears 
that the main resistance to the gravitational pull of the matter fields comes from the magnetic 
tension. The latter reflects the elasticity of the field lines and their tendency to remain straight. 
This tendency is manifested through a variety of tension stresses that try to restore equilibrium 
every time the magnetic forcelines are distorted. 

General Relativity describes the macroscopic behaviour of the gravitational field. In the 
geometrical framework of the theory, gravity is represented by the curvature of the spacetime, 
while the interplay between matter and geometry is monitored by Einstein's equations. However, 
for energy sources of vector nature, like the electromagnetic field for example, an additional 
set of equations is of almost equal importance. These are the Ricci identities. Our analysis 
shows that, within framework of General Relativity, one of the agents causing magnetic-line 
deformations is the curvature of the space. The mathematics of the interaction are monitored 
by means of the Ricci identities [17J. As the collapse proceeds and the departure from the 
Euclidean flatness increases, the aforementioned deformations grow as well. As a result, the 
restoring tension stresses also get stronger. This opens the theoretical possibility that, when the 
curvature distortion crosses a certain threshold, the resulting tension stress may become strong 
enough to stop the contraction and, perhaps, turn it into expansion [17J. Just like a spring 
bursts forward after been compressed back to its limit. We revisit this scenario in § 14.31 14.41 
and show that some degree of anisotropy would probably be necessary, if the magnetic tension 
is to bring the collapse to a halt. The required amount of anisotropy is found to depend on the 
.B-field's relative strength. We attempt to quantify these statements by providing the necessary 
conditions for the tension stresses to counterbalance the gravitational pull of the matter in 
Raychaudhuri's equation. 
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2 The 1+3 covariant approach 



The covariant approach to general relativity utilizes an 1+3 'threading' of the spacetime, as 
opposed to the 3+1 'slicing' of the ADM decomposition (see [27 1 for a recent review). This is 
achieved by introducing the timelike 4- velocity u a = dx a /dr, normalised so that u a u a = — 1. 
The it a -field defines the worldlines of the fundamental observers (with r being the associated 
proper time), while the tensor h ab = g ab + u a u b (with g ab representing the spacetime metric) 
projects orthogonal to u a and into the observers' 3-dimensional, instantaneous rest-space. On 
using u a and h ab , every variable, operator and equation spits relative to the ii a -field into their 
irreducible timelike and spacelike components. 

2.1 The gravitational field 

Within the geometrical interpretation of general relativity, the gravitational field is described by 
means of the Riemann curvature tensor. The latter splits into a local (Ricci) and a long-range 
(Weyl) part according to 

Rabcd = Cubed + ~ (dacRbd + QbdRac ~ QbcRad ~ dadRbc) ~ {QacQbd ~ 9ad9bc) , (1) 

2 o 

where C abc d is the (conformal curvature) Weyl tensor and R ab is the Ricci tensor (with R = R a a 
being the associated Ricci scalar). The Riemann tensor also obeys the symmetries R a bcd = 

R[ab][cd] = Rcdab- 

The Ricci field is determined by the local matter distribution and the interplay between 
geometry and matter is governed by the Einstein field equations. In the absence of a cosmological 
constant the latter read 

Rab ~ ^ R 9ab = KT ab , (2) 

where k = 8ttG and T a b is the (symmetric) energy-momentum tensor of the matter. This is 
determined by the distribution of the local matter and satisfies the conservation law V 6 T a 6 = 0. 

While the Ricci tensor monitors the local gravitational field, the Weyl tensor describes its 
long-range counterpart. In other words, the conformal curvature tensor represents the part of 
the spacetime curvature that is not determined by the local matter, but depends on the matter 
distribution at other points. This relates the Weyl field to tidal forces and gravitational waves. 
Further insight is obtained when C a b c d is decomposed into its irreducible constituents. Given a 
4- velocity congruence u a , we have 

C ab cd = 4 (u [a u^ + h [a ^ E b f + 2e abe u^H d ^ + 2u [a H b]e e cde , (3) 

with E ab = C acb dU c u d and H ab = e a cd C c d be u e /2 representing the electric and magnetic parts of 
the Weyl tensor (relative to the u a -frame). Note that the magnetic component has no Newtonian 
analogue and for this reason is more closely (though not explicitly) associated with gravitational 
radiation. Also, both parts of the Weyl field are necessary if gravitational waves are to exist. 
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2.2 Kinematics 



An advantage of 1+3 decomposition is its mathematical compactness and physical transparency, 
since every single one of the irreducible components involved has a very clear physical or geo- 
metrical interpretation. For example, the covariant derivative of the 4-velocity field splits as 

V b u a = D b u a - A a u b = - Qh ab + a ab + u ab - A a u b , (4) 

where D a = h a b V b is the 3-D covariant derivative operator. The right-hand side of the above 
provides the irreducible kinematic variables that describe the observers' motion. These are 
the volume scalar, = V a u a = D a u a , the shear and the vorticity tensors, a ab = D( b u a j and 
u ab = D[ b u a ] respectively, and the 4-acceleration vector, A a = ii a u The latter represents non- 
gravitational forces and therefore vanishes when matter moves along geodesic worldlines. Also, 
a positive/negative volume scalar implies expansion/contraction for the fluid element, the shear 
modifies its shape and the vorticity changes its orientation. Note that, starting from the vorticity 
tensor, we may also define the vorticity vector by means of oj a = e abc oj bc /2, with e abc representing 
the (totally antisymmetric) 3-D Levi-Civita tensor. 

The motion of the fundamental observers is described by a set of three propagation equations 
supplemented by an equal number of constrains. These follow after applying the Ricci identities, 

2V[ V 6 ]« C = R a bcdu d , (5) 

to the associated 4-velocity field. The timelike component of the above leads successively to 
Raychaudhuri's formula, 

® = -\® 2 -\<P + 3p) -2 {a 2 - lo 2 ) + D% + A a A a , (6) 

to the shear-propagation equation 

2 1 

°(ab) = -g ®Vab ~ <?c(a° C b) ~ ^(a^b) + ^>(a A b) + A (a A b) ~ E ab + - KTT ab (7) 

and to the evolution equation of the vorticity 

2 1 

^(a) = - g - - curL4 a + a ab u h , (8) 

where a 2 = a ab a ab /2 and to 2 = u ab io ab /2 = uj a u a . Expression ([6]) dictates the average kinematics 
of the fluid element, while (JT]) and © monitor the evolution of kinematically induced anisotropies 
and the rotational behaviour of the fluid respectively. For our purposes, the key equation is 
Raychaudhuri's formula, which governs the rate of the volume contraction/expansion of a self- 
gravitating medium. Negative terms in the right-nahd side of © assist the collapse (or inhibit 
the expansion), while positive ones resist contraction (or lead to expansion). 

1 Round brackets denote symmetrisation and square ones antisymmetrisation. Angular brackets are used to 
describe the symmetric and trace-free part of orthogonally projected second rank tensors and vectors (e.g. a a t = 
D(bWa) = D(i,« a ) — (D c u c /3)h a b and E( a ) = h a b E a ). Also, overdots indicate covariant differentiation with respect 
to proper time, namely along the fundamental worldlines (e.g. A a = ii a = u b \7bU a )- 
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The spacelike component of Eq. ([5]) also splits in an analogous way. The result is expressed 
ELS cl set of three (kinematic) constraints, given by 

2 

DVafc = - D a G + curlw a + 2e abc A b uj c - nq a , (9) 
D a u a = A a uj a (10) 

and 

H ab = curkr afo + D( a w 6 ) + 2A( a u b ) , (11) 

respectively. Together, relations ([6])- (jlip provide a fully relativistic kinematic description of a 
general imperfect fluid. 

2.3 Matter fields 

The energy- momentum tensor of the matter also decomposes into its irreducible parts. When 
dealing with an imperfect fluid, for example, we have 

Tab = PU a U b + phab + 2<?( a U b ) + 1T ab . (12) 

Here, p = T ab u a u h is the energy density, p = T ab h ab /3 is the isotropic pressure, q a = —h a b T bc u c 
represents the energy flux and 7r ab = h( a c hfy d T C( i the anisotropic pressure, all measured in the 
fundamental frame. In the case of a perfect fluid, q a and ir ab vanish identically and the above 
reduces to 

Tab = PU a U b + phab ■ (13) 

with p and p related through the equation of state. For barotropic media, the latter takes the 
simple p = p{p) form. 

The nature of the matter fields also determines the conservation laws. Assuming a general 
imperfect fluid, the timelike component of the (twice contracted) Bianchi identities provides the 
energy conservation law 

p = -@(p +p)- B a q a - 2A a q a - a ab TT ab • (14) 
Its spacelike counterpart, on the other hand, leads to the conservation of the momentum 

4 b b b 

(p + p)A a = ~B a p ~ q(a) ~ g @<?a ~ {<Jab + ^ab) Q ~B 7T ab - 7T a bA . (15) 

As mentioned above, when dealing with perfect fluids, we may set q a = = 7r ab . Then, expres- 
sions (114j) and (|15h simplify considerably. 

3 Electromagnetic fields 

The electromagnetic field is a source of energy density and, as such, it contributes to the right- 
hand side of the Einstein equations (see (J2])). In addition, the Maxwell field has its own evolution 
and constraint formulae, while it also obeys a separate set of conservation laws. 
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3.1 Fluid description 

The electromagnetic field is invariantly described by the antisymmetric Faraday tensor (F^), 
which relative to our fundamental observers, splits as 

F ab = 2u [a E b] + e abc B c (16) 

where E a = F ab u b and B a = e abc F bc /2 are the associated electric and magnetic fields respectively. 
The above decomposition also facilitates a fluid description of the Maxwell field, analogous to 
that seen in (|12p . To be precise, given a 4- velocity field, the electromagnetic energy-momentum 
tensor splits as 

T ab m) = \i E2 + B 2 ) u « u b + ~ {E 2 + B 2 ) h ab + 2Q {a u b) + U ab , (17) 

where E 2 = E a E a and B 2 = B a B a are the square magnitudes of the individual fieldsH Also, 
Qa = £abcE b B c is the electromagnetic Poynting vector and II a b = —E^E^ — B( a B b y Comparing 
the above to Eq. (|12p . confirms that the Maxwell field corresponds to an imperfect fluid with 
energy density p( ern > = (E 2 +B 2 )/2, isotropic pressure p( em ) = (E 2 +B 2 )/6, energy flux q e a m = a 
and anisotropic stresses given by H a b- Finally, decomposition (|17j) ensures that ri em)a =o, 111 
agreement with the trace-free nature of the electromagnetic radiation. 



3.2 Maxwell's equation 

Invariantly, Maxwell's equations are given by the set V[ a F bc j = and V b F ab = J a , with J a rep- 
resenting the electric 4-current. On using decomposition (|16p . the above expressions decompose 
into a set of two propagation equations 

E {a) = --QE a + (a ab + e abc oj c )E b + £abc A b B c + cuvlB a -Za, (18) 



B {a) = --QB a + (a ab + e abc co c )B b -e abc A b E c -cmlE a , (19) 
supplemented by the constraints 

D a E a + 2co a B a = n (20) 

and 

B a B a - 2u a E a = . (21) 

Note that curll? a = e abc D b B c by definition, with an analogous expression for curlE^. Also, Za 
represents the orthogonally projected (spatial) electric current and \x is the electric charge. The 
two are related to the 4-current by means of 3a = h a b J b and \i = —J a u a respectively. As a 
result, we may write 

J a = nu a + 3 a . (22) 

Relations (|18p -([2TT) are the covariant analogues of Ampere's law, Faraday's law, Coulomb's 
law and Gauss' law respectively. A characteristic feature of these formulae is that they incor- 
porate relative-motion effects, carried by the kinematic terms in their right-hand side. The first 



2 We use Heaviside-Lorentz units for the electromagnetic field throughout this article. 
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two terms in the right-hand side of the Ampere law, in particular, represent the magnetic field 
induced by the relative motion of its electric counterpart. Analogous terms are also seen in Fara- 
day's law. Similarly, kinematical effects introduce an effective electric charge and an effective 
magnetic charge in Eqs. (|20p and (|2ip respectively. As a result, the magnetic vector is generally 
not a solenoidal. Finally, the acceleration terms in expressions (|18|) and (|19|) reflect the fact that 
spacetime is treated as a single, unified entity. 

3.3 Conservation laws 

The electromagnetic stress-energy tensor is separately conserved (i.e. \/ b T^ m ^ =0). In the 
pesence of charged matter, however, one needs to substitute the total energy-momentum tensor 
into the Bianchi identities and then decompose the resulting relation. Then, using the full 
Einstein-Maxwell system, the timelike part of the aforementioned splitting reads 

p = -G(p + p) - B a q a - 2A a q a - a ab 7r ab + E a Z a , (23) 

and provides the energy conservation law. The spacelike component, on the other hand, leads 
to the conservation of the momentum 

(p + p)A a = -B a p - q {a) - - @q a - (a ab + u ab ) q b - Tl\ ab - Tr ab A b + pE a + e abc Z b B c . (24) 

Note that q a = = 7r ab for perfect fluids, in which case the right-hand sides of ([23]) and ([23]) 
simplify considerably. Also, as we shall see below, the overall effect of the electromagnetic terms 
in the above depends on the electric properties of the matter. 

Finally, when dealing with charged media, there is a an additional conservation law. This 
is related to the 4-current and is given by the constraint V a J a = 0. The latter leads to the 
conservation law of the electric charge 

A = -®M " D a X - A a Z a , (25) 
expressed relative to a family of fundamental observers. 

3.4 Ohm's law 

The evolution and the effects of the Maxwell field depend on the electrical properties of the 
charged medium, which are described by Ohm's law. At the resistive MHD limit and in the 
frame of the fluid, the latter takes the covariant form [28\ [29] 

J a = fJ-u a + sE a , (26) 

with ^ representing the (scalar) conductivity of the medium. Substituting decomposition (I22D 
into the left-hand side of ([26]) and then projecting orthogonal to u a , gives 

Za = sE a . (27) 

This, alternative form, of Ohm's law relates directly the 3-current to the S-field through the 
electrical properties of the host medium. 
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The two limiting cases are those of very poor and very high electrical conductivity, corre- 
sponding to ? — > and ? — > oo respectively. Following Eq. (|27p . the former has zero electric 
currents and the latter is characterised by a vanishing electric field. For any other intermediate 
situation, one needs a model for the conductivity of the matter. 

4 Charged collapse 

Over the years, the gravitational contraction of charged matter (usually dust) in the presence 
of an electromagnetic (usually a purely electric) field has been considered by several authors 
from a variety of perspectives (see [l]-[20j for a representative though incomplete list). A com- 
mon feature in a number of these studies, was the conclusion that charged dust could, under 
'favourable' circumstances, escape the ultimate collapse to a central singularity. 

4.1 Poorly conductive dust 

Some of the aforementioned approaches are static, others are dynamical and, with very few 
exceptions, all are metric-based. Here, we will take a dynamical look into the question of 
charged collapse. We will also approach the problem from the 1+3-covariant point of view, 
instead of adopting metric-based techniques. This puts the Raychaudhuri equation (see relation 
© in § 12. 2p at the centre of our analysis. 

Raychaudhuri's formula has played a fundamental role in gravitational-collapse studies and 
has been at the centre of all the singularity theorems. Expression ([6]) is the general one and can 
be adapted to any physical environment. When dealing with charged matter in the presence of 
an electromagnetic field, one needs to take into account the fact that the Maxwell field adds to 
the gravitational mass of the system and also incorporate the contribution of the Lorentz force 
to the momentum-conservation law. Assuming matter in the form of a pressureless perfect fluid, 
we may set p, q a and ir a b to zero. Then, at the low conductivity limit, Za vanishes as well and 
the momentum-conservation law (see expression (I24D in § 13. 3p reduces to 

A a = £E a , (28) 

where £ = p/p measures the specific charge density of the fluidH Note that this ratio is time 
invariant in our case. Indeed, for a pressureless perfect fluid and for poor electrical conductivity, 
both p and p obey identical evolution laws (see Eqs. (|23l) and (i25j) in § I3.3j) . We also underline 
the absence of any direct magnetic effects in the above relation, which results from the vanishing 
of the 3-currents. 

In order to analyse the implications of a nonzero net charge for the average kinematics 
(i.e. the volume expansion/contraction) of a self-gravitating fluid, we need to substitute (|28p 
into the right-hand side of Eq. ([6]). Then, on using Coulomb's law (see expression (I20p in § 13. 2p . 
the latter reads 

9 = - l -Q^- l -n{p + E 2 + B 2 )-2{a 2 -u 2 )+e{p + E 2 ) 

+E a T) a i - 2& a B a . (29) 

3 To the best of our knowledge, all the available studies of charged collapse also ignore the presence of the 
electric currents. Seen it from our perspective, this corresponds to choosing poorly conductive environments. 
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In contrast to Eq. (|28p . the above contains direct magnetic effects. These come from the B- 
field's contribution to the total energy density of the system and also to the effective charge (see 
relation (|2Uj) in § I3.2p . For our purposes, however, the key quantity is the (positive definite) 
fourth term in the right-hand side of Eq. (|29p . This term reflects the presence of a net electric 
charge and monitors the effect of the repulsive Coulomb forces on our self-gravitating system. 

4.2 Coulomb repulsion vs gravitational attraction 

Assuming that spherically symmetry and spatial homogeneity are largely preserved, we may 
ignore the shear and the vorticity, as well as variations in the spatial distribution of the charge 
density. In that case Raychaudhuri's equation simplifies to 

® + \& 2 = -\k{ P + E 2 + B 2 ) + ? {p + E 2 ) . (30) 

The above describes the average volume contraction/expansion of an almost spherically symmet- 
ric, charged perfect fluid of low electrical conductivity!!! Put another way, Eq. ([30]) monitors the 
average separation between the (timelike) worldlines of two neighbouring observers that move 
along with the matter. We remind the reader that, in the presence of a nonzero net charge, 
these worldlines are no longer geodesies. 

In the right-hand side of (|30p one can clearly see the competing action of the gravitational and 
the Coulomb forces and how these affect the kinematics of the charged fluid. The first term, in 
particular, shows how the Maxwell field increases the overall pull of gravity by adding to the total 
energy density of the system. Alternatively, one might say that the electric charge contributes to 
the total gravitational mass, an effect occasionally referred to as 'charge regeneration' pQ. The 
second term, on the other hand, demonstrates how the repulsive Coulomb forces act against 
gravity. Note that this stress is partly Newtonian and partly relativistic. The Newtonian 
component is £ 2 p and derives from the 3-divergence of the 4-acceleration in Eq. ([6]), which 
also appears in non-relativistic treatments |30|. 13 lj . The relativistic effects, on the other hand, 
are encoded in the £ 2 E 2 stress. The latter comes from the square of the 4-acceleration in the 
right-hand side of the Raychaudhuri equation and has no Newtonian counterpart [301 13T) . All 
these mean that the electric field contributes equally to the total gravitational field and to the 
Coulomb forces. In contrast, the magnetic field only adds to the overall pull of gravity. This 
difference reflects the absence of magnetic monopolesH 

It is a very well known result that, as long as the right-hand side of Eq. (I30p remains non- 
positive, an initially converging pair of worldlines will focus within finite proper time. This 
implies that, within our approximation scheme, the Coulomb forces can stop these worldlines 
from focusing (and therefore prevent caustic formation) if 

e(p + E 2 )> l -K{p + E 2 + B 2 ) . (31) 

4 In the case of exact spherical symmetry the magnetic field vanishes in the frame of the fluid, while only the 
radial component of its electric counterpart survives (e.g. see [16U18| ). 

5 If magnetic monopoles were also included into the scheme, their presence would have led to a Coulomb-like 
repulsive stress in the right-hand side of Eq. (|30|) . analogous to the one produced by the electric charges. 
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When the above condition holds, the worldlines of a contracting, nearly spherically symmetric, 
poorly conductive, pressureless, charged fluid will not focus. Moreover, in the absence of the 
magnetic field, as it happens in the case of exact spherical symmetry for example, the above 
reduces to 

t>]fY K = ]fl M ™> (32) 

having set k = Mp 2 (with Mpi representing the Planck mass). Note that the last condition has 
also been discussed, as a requirement for avoiding both the initial and shell-crossing singular- 
ities during the spherically symmetric collapse of charged dust, in the metric based approach 
of [16] . Analogous conclusions, relating the charge-to-matter ratio with the stability (static or 
dynamical) of a spherically symmetric charged medium against gravitational contraction, were 
also reached in the earlier treatments of [H [9] . 

The final form of ([31]) also depends on the strength of the electromagnetic components 
relative to the matter content. Assuming that E 2 , B 2 -C p, we arrive again at £ > a/ k/2 = 
Mpi /\/2. When E 2 ~ B 2 » p, on the other hand, condition ([3T]) recasts into £ > yJH = Mp\. 
Finally, for E 2 ~ B 2 ~ p, expression (j3~Tj) takes the form £ > y / 3At/4 = \/3/ 4 Mpi . In all these 
cases the repulsive Coulomb forces will be strong enough to counterbalance the gravitational 
field and prevent the (non-geodesic) worldlines of the charged particles from focusing. 

Before closing this section, we should also briefly discuss the subtle role played by possible 
inhomogeneities in the spatial distribution of the electric charge. These are incorporated in the 
second-last term of Eq. ([29]) and their effects have so far been bypassed. When included into 
expression ([30]) . the latter reads 

9 + i G 2 = ~k {p + E 2 + B 2 ) + £ 2 [p + E 2 ) + E a D a £ , (33) 

where E a D a ^ ^ in general. Depending on its sign, the inhomogeneous term in the right-hand 
of the above will either strengthen or weaken conditions (|3ip and (|32p . In particular, when 

E a Vai = £ f - E a B a p - - E a B a p) > , (34) 
\H P J 

the spatial gradients of the specific charge will assist the Coulomb repulsion and therefore hamper 
the collapse further. In the opposite case, the situation is reversed and the inhomogeneous term 
in the right-hand side of (I33p enhances the overall gravitational pull. This versatility in the 
effect of the charge gradients has been noted in earlier studies as well [8] [19] . The reader is 
particularly referred to [19] for an alternative, metric-based, discussion of the matter. 

4.3 Highly conductive dust 

The assumption of a charged-matter distribution in the absence of magnetic fields and electric 
currents is probably a highly idealised one. Put another way, poorly conductive environments 
offer a rather unlikely description of stellar interiors. Strong magnetic fields and high electrical 
conductivity sound much more likely. 

When dealing with highly conductive matter, Ohm's law (see Eq. ([27]) ) guarantees that the 
electric fields vanish in the frame of the fluid. At the same time, the finite currents keep the 
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magnetic field lines frozen in with the matter [32\ [33] . This is the familiar ideal-MHD approx- 
imation, which has long been used in studies of magnetised plasmas in both astrophysics and 
cosmology. Here, we will also assume overall charge neutrality. In such environments, Maxwell's 
formulae (see relations (|18p -(|2i p ) reduce to one propagation equation and three constrains. More 
specifically, Ampere's law leads to the magnetic induction equation 

B {a ) = -3 &B a + (a ab + e abc oj c )B b , (35) 
while Faraday's, Coulomb's and Gauss' laws reduce to the constrains 

e abc A b B c + cm\B a =Za, u a B a = and B a B a = 0, (36) 

respectively. Note that Eq. (|35p guarantees that the magnetic lines always connect the same 
particles (i.e. that the -B-field is frozen-in with the matter). 

The kinematics of a highly conductive pressureless fluid are still governed by the set ©- 
(fTTI) . once the latter is adapted to the specifics of the MHD limit. In particular, Raychaudhuri's 
equation reads 

6 = -i 9 2 - i k (p + B 2 ) - 2 (a 2 - oj 2 ) + B a A a + A a A a , (37) 

with the 4-acceleration given by the momentum-conservation law, which has now changed to 

(p + B 2 ) A a = -e abc B b cmlB c = - ~ D a B 2 + B b B b B a , (38) 

The above set monitors the average contraction/expansion of a self-gravitating, highly conduc- 
tive, pressureless fluid in a magnetised environment. At this point, we would like to draw the 
readers attention to the second equality in Eq. ()38|) . where the Lorentz force has been split into 
its pressure and tension components. The former comes from the positive pressure that the 
-B-field exerts orthogonal to its own direction and demonstrates the tendency of the magnetic 
forcelines to push each other apart. The latter reflects the negative pressure, namely the tension, 
that the field exerts along its own direction. This stress manifests a unique feature of magnetic 
fields: the elasticity of their forcelines and their tendency to remain 'straight' \32\ I33j. 

4.4 Magnetic tension effects 

Because of their elastic properties, the field lines react to any attempt that distorts them from 
equilibrium. In standard Newtonian MHD, the agents causing such deformations are ordinary 
charged particles. The field's reaction then triggers a variety of tension stresses, which in turn 
lead to a range of generally well studied effects. When General Relativity is involved, however, 
there is an additional agent that can affect the distribution and the topology of the magnetic 
forcelines. This new agent is the curvature of the spacetime itself. 

The vector nature of magnetic fields means that they 'feel' the curvature of the space via 
the Ricci identities, as well as through Einsten's equations. The former are a purely geometrical 
set of relations that have played a fundamental role in the mathematical formulation of general 
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relativity, since they essentially provide the definition of spacetime curvature. Applied to the 
magnetic vector, the Ricci identities read 

2V [a V b] B c = R abcd B d , (39) 

with R a bcd representing the Riemann tensor of the spacetime (see Eq. (H|) in § 12 . X [) . Starting 
from the above, one can obtain the so-called 3-Ricci identity 

2D[ a Db]-Bc = lZ a bcd,B d — 2co ab B^ , (40) 

which now involves TZ abc( i, namely the Riemann tensor of the 3-dimensional hypersurfaces or- 
thogonal to our fundamental observers. Note the vorticity term in the right-hand side of the 
above, which reflects the fact that in a rotating spacetime the observers' worldlines are not hy- 
persurface orthogonal. The Ricci identities demonstrate, mathematically, how distortions in the 
geometry of the host spacetime translate into magnetic- line deformations. These in turn lead to 
a range of tension stresses that are purely geometrical in origin. Next, we will concentrate on 
the implications of such tension stresses for the gravitational collapse of a magnetised fluid. 

For our purposes, there is no real loss of generality in assuming (almost) zero rotation and 
(nearly) spatially homogeneous distributions for both the matter and the magnetic energy densi- 
ties. The latter assumption, in particular, allows us ignore the pressure part of the Lorentz force 
(see Eq. (I38p ) and focus on its tension component. In that case, the momentum-conservation 
law reduces to (p + B 2 )A a = B b D b B a . Substituted this expression into the right-hand side of 
Eq. (|371) and using relations (|36b ) , (fiOj) , gives 

6 + i G 2 = -R ab u a u b + c 2 a TZ ab n a n b - 2 (a 2 - a\) , (41) 

where R ab and 1Z ab are the 4-dimensional and the 3-dimensional Ricci tensors respectively c 2 = 
B 2 /(p + B 2 ) is the square of the Alfven speed and n a is the unit vector along the direction of 
the magnetic field. Also note that 2a 2 B = (D^ a B b ^) 2 /(p + B 2 ) and we have dropped the last 
term of Eq. (|37[) . That quantity, namely the square of the 4-acceleration - A a A a , always resists 
gravitational contraction and its absence will not weaken our argument (see below). 

The first term in the right-hand side of Eq. (|4ip provides the total gravitational mass of the 
system, with R ab u a u b = n{p + B 2 )/2 > in our case. Typically, this is the main quantity that 
drives the contraction, with the shear adding to the overall pull of gravity. The magnetic terms 
in (I41I) are due to the field's tension and are triggered by the deformation of its forcelines. Both 
terms represent restoring stresses, which act against the agent that caused the line distortion in 
the first place. The magneto-shear, for example, is triggered by shearing effects. This explains 
why gb acts against shear proper in Eq. (|41f) . Note that this term is non-relativistic in origin, 
since it also appears in purely Newtonian treatments [31jl^l 

6 When rotation is also included into the system, it adds a tension stress to the right-hand of I41 1) that may 
be seen as the magnetic analogue of vorticity proper 17 . This stress is triggered by the twisting of the field 
lines, which trapped within the rotating fluid. In analogy with the magneto-shear term in Eq. (|41[) . the magnetic 
vorticity stress acts against its kinematic counterpart, an effect that is sometimes referred to as 'magnetic braking'. 
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4.5 Magnetic tension vs gravitational attraction 

Of all the quantities in the right-hand side of (|41|) . the most unusual and intriguing is the second. 
This also represents a tension stress, but it is purely relativistic in nature. The difference with 
the other tension stresses is that, now, the agent responsible for the deformation of the magnetic 
field lines is the geometry of the 3-space. In particular, spatial curvature distortions affect the 
topology of the magnetic forcelines, giving rise to the mageto-curvature tension stress seen in 
Eq. (|4ip . As mentioned earlier, the whole interaction is mathematically monitored via the Ricci 
identities. Not surprisingly, the magnitude of the tension stress depends on the relative strength 
of the -B-field and on the amount of the line deformation. The latter is determined by the sum 
TZ ab n a n b , which represents the 3-curvature distortion along the direction of the field lines. What 
is surprising, is that the effect of the magneto-geometrical stress in (I4ip depends on the sign 
of the curvature deformation. Thus, the aforementioned stress will resist the collapse when the 
scalar lZ ab n a n b is positive, something that is generally expected in all cases of realistic stellar 
collapse. Moreover, the strength of the tension stress increases with growing curvature distortion, 
just like the resisting power of a spring grows when it is compressed back to its limit. All these 
mean that, at least in principle, it is conceivable that the aforementioned stress can grow strong 
enough to halt the contraction and therefore prevent the particle worldlines form focussingHI 
Whether this happens or, depends on a number of parameters. Perhaps the decisive factor is 
the balance between the two geometrical terms in the right-hand side of Eq. (I4ip . Whether, 
in particular, the tension stress can outbalance the gravitational pull of the matter. In what 
follows we will take a closer look into this possibility. 

Suppose that at some time during the collapse of a magnetised fluid the magneto-geometrical 
tension stress has grown strong enough, so that 

c 2 a K ab n a n b > RabU a u b . (42) 

Then, substituting the Gauss-Codacci equation (e.g. see expression (1.3.39) in [57] - with zero 
rotation in our case) into the left-hand side of the above, we obtain 

| Up ~ ~ © 2 ) + {E ab - i ®<J ab + <J ca a c b ^j n a n b > ~ k (p + B 2 ) . (43) 

The first of the two parentheses on the left represents the isotropic part of the tension stress 
and the second the anisotropic. In particular, E ab represents the electric component of the Weyl 
tensor, which is known to describe the free gravitational field. We now introduce the auxiliary 
variables f3 = B 2 / p and Cl p = Snp/Q 2 . The former measures the relative strength of the -B-field 
and the latter corresponds to the familiar density parameter typically associated with the FRW 
universes. Then, condition ()43p takes the form 

Ka b n a n b > 1 k P (l + ~ p + | /T 1 + 2UjA > , (44) 

where d\ ab = E ab — (Q/3)a ab + o ca o c b . The fact that the scalar 9\ ab n a n b is nonzero implies that 
the collapse must be anisotropic, if the tension stress is to outbalance the gravitational pull of 



7 Recall that, in the presence of a non-zero Lorentz force, the worldlines in question are no longer geodesies. 
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the matter. Calculating the required amount of anisotropy is a rather involved process, as it 
depends on a number of (essentially free) parameters. Nevertheless, some general estimates are 
still possible. 

The first, qualitative, result comes from the dependence of < 3\ a b'n a n b on the electric Weyl 
tensor and also on the shear. This dependence suggests that the tension stresses are more 
likely to outbalance the pull of gravity in situations where the Weyl part of the gravitational 
field dominates over its Ricci counterpart. More quantitative estimates can be obtained by 
assuming that Q, p S> 1 during the advanced stages of the contraction, which sounds plausible 
for most realistic collapse scenarios. Then, following (|44|) . the required amount of anisotropy 
depends entirely on the relative strength of the -B-field (i.e. on the /3-parameter) . Assuming that 
the two sources are equally strong (i.e. setting /3 ~ 1), we find that condition (j44l) is satisfied if 
< y\ a i J n a "n b > p. This changes to 9\ a i,n a n b > B 2 when the magnetic field dominates (i.e. for ^ > 1). 
In the opposite case, namely when {3 <C 1, expression ()44p leads to 9\ a bn a n b 3> p. Overall, 
when the magnetic field is equally strong or stronger than the matter, the tension stresses can 
outbalance the gravitational pull of the matter fields provided the collapse is sufficiently (though 
not necessarily highly) anisotropic^ This changes when the -B-field is weak. Then, the magneto- 
geometrical tension stress will remain subdominant unless the collapse is highly anisotropic. 
Note, however, that the anisotropy requirements may weaken (perhaps significantly), if the 
effect of the repulsive j4 a ^4 a -stress in the right-hand side of the Raychaudhuri equation (see 
expression (|37l) ) are also accounted for. 



5 Discussion 



The gravitational contraction of charged matter and the implications of electromagnetism for 
the fate of the collapse have been studied by many authors, from a variety of perspectives and 
by means of different methods. A number of these studies raised the possibility of a 'bounce', 
which sometimes was caused by the presence of a net charge, while in others it was triggered 
by the presence of the magnetic field. The available approaches are static as well as dynamical 
and in almost all cases they employ metric-based techniques. 

In this work we have adopted covariant methods to pursue a dynamical study of charged 
collapse. This has put the Raychaudhuri equation at the centre of our analysis. We have 
considered two cases that correspond to two opposite environments. The first assumes poorly 
conductive matter and looks at the effects of net charges, and on those of the resulting repulsive 
Coulomb forces, upon gravitational contraction. The second considers magnetised environments 
of high electrical conductivity and focuses on the role of the magnetic tension. In the first 
case we found that, to a very large extent, the overall outcome of the collapse depends on the 
value of the specific charge (i.e. on the ratio between the charge and the energy densities) of 
the medium. When the charge density dominates, the Coulomb forces become stronger that the 
gravitational pull of the matter and the contraction stops. In the opposite case, the collapse 
proceeds essentially unimpeded by the Coulomb stresses. Our results agree both qualitatively 

8 Even when the magnetic field dominates (i.e. when B 2 3> p), the total fluid satisfies the three energy conditions 
(i.e. weak, strong and dominant). One can easily verify this by recalling that ps = B 2 /2, while pb = —B 2 /2 
along the direction of the field lines (i.e. where the magnetic tension acts) and pb = B 2 /2 orthogonal to them. 
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and quantitatively with those obtained by earlier, primarily metric-based, studies. 

Highly conductive environments probably offer a better approximation of realistic gravita- 
tional contraction. This alternative corresponds to the familiar ideal-MHD limit, where there 
are no electric fields and the currents keep the magnetic lines frozen- in with the fluid. In the 
absence of net charges, we found that the main resistance against gravitational collapse is prob- 
ably the one coming from the magnetic tension. The latter reflects the elasticity of the field lines 
and their tendency to react against anything that distorts them from equilibrium. Assuming 
that General Relativity is the correct theory of gravity, one of the agents leading to such line- 
deformations is the curvature of the space. This means that magnetic forcelines trapped within 
a collapsing medium are being increasingly deformed by the ever growing gravitational field. 
These deformations inevitably lead to magneto-geometrical tension stresses that resist further 
contraction. A close analogy is that of a spring's tension, which grows stronger as the latter 
is compressed back to its limit. It is therefore conceivable that, under certain circumstances, 
the aforementioned elastic stresses may become strong enough to outbalance the pull of gravity 
and bring the collapse to a halt. Our analysis shows that, if something like that is to happen, 
it would probably need some degree of anisotropy. Although the exact amount of the required 
anisotropy is rather difficult to estimate at the moment, Weyl-curvature dominated spacetimes 
may offer a suitable host to test this hypothesis. 
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